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<£ ; Abstract 

<N 

Uncertainties in the theoretical prediction for the inclusive B — > X s 7 decay 

rate are discussed. We emphasize that there is no operator product expansion 
(NJ I for this process. Nonetheless, some nonperturbative effects involving a virtual 

cc loop are calculable using the operator product expansion. They give a 
pg ■ contribution to the decay rate that involves the B meson matrix element of 

an infinite tower of operators. The higher dimension operators give effects 

that are only suppressed by powers of ttt.;,Aqcd/wIc ~ 0.6, but come with 

small coefficients. 



The inclusive B — > Xg'j decay has received considerable attention in recent years [1-15], 
since it is sensitive to physics beyond the standard model || (X s denotes a final state with 
strangeness —1 and charm 0). The photon spectrum also carries information on nonper- 
turbative physics that can help us better understand other B decays |13| , p^] . The recent 
CLEO measurement |1J excludes large deviations from the standard model. Therefore, it is 
important to know the standard model predictions as precisely as possible. 

Since the b quark is heavy compared to the QCD scale, one would hope that the inclusive 



B —* X s 7 decay rate can be calculated in a systematic QCD-based expansion [|16 |. The 
dominant contribution to the decay rate comes from the matrix element of the electromag- 
netic penguin operator (usually denoted by 7 ). In the m b — > oo limit, it is given by the 
free quark decay result. The leading nonperturbative corrections to this contribution are 
suppressed by (Aq CD /mb) 2 . Provided the photon energy is not restricted to be too close to 



its maximal (i.e., end-point) value, they are quite small, around —3% [|12 ]. With the recent 



completion of the full next-to- leading order perturbative calculation fLlfl , it is usually argued 
that theoretical uncertainties in the prediction for the inclusive B — > X s ^ decay rate are 
not larger than 10%. 

The effective weak interaction Hamiltonian at a scale /i (of order m^) is given by 

AC 8 

tfe ff = --£ v* v* E c*(aO Oi(ji) ■ (i) 

In the conventional notation, 2 = {sLal^L^) (cl@ I^clo), Ox om y differs from 2 in the 
way color indices are contracted, O3 — Oq are four-quark operators involving all flavors below 
the scale /x, O-j = (e/167T 2 ) m& s^ cr^F^ 6jj, and Og is obtained from O-j by replacing eF F 
by QsG^y. Using perturbative QCD to evaluate the B —> Xgj decay rate, in the leading 
logarithmic approximation the matrix element of Cy(/x) 0-?(n) dominates for large enough 
photon energies. 

A systematic computation of the Oj contribution to the inclusive B — > X s 7 decay rate 
involves performing an operator product expansion (OPE) for the time ordered product 

T 77 = ^- J d*xe~^ (B(v)\T{Of(x)O^(0)}\B(v)) g, u , (2) 




FIG. 1. The integration contour C in the complex v-q plane. The cuts extend to Hev-q — ► ±00. 

to all orders in the strong interaction. Here O7 = (i e/8n 2 ) mi, sl u^qx 6r. At fixed q 2 = 0, 
this time ordered product has a cut in the complex v-q plane along v-q < mb/2 corresponding 
to final hadronic states X s , and another cut along v ■ q > 3m&/ '2 corresponding to final 
hadronic states X^. The contribution of the magnetic moment operator O7 to the B — > X s 7 
decay rate is given by the discontinuity across the cut in the region < v ■ q < m&/2, 

dr AG 2 F \V t * s V tb \ 2 C 2 
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Since the cuts are well-separated, one can compute this contribution to the B — ► X s 7 decay 
rate assuming local duality at the scale m^. (The integration of T 77 over the contour C in 
Fig. 1 pinches the physical cut at v • q — 0, but at that point the hadronic final states have 
invariant mass mx s = ttlb-) 

At leading order in the OPE, the dimension-three operator 67^6 occurs. Its matrix 
element gives a calculable contribution to the photon energy spectrum proportional to 5(E 1 — 
171^/2). Higher dimension operators give terms proportional to derivatives of this delta 
function, and the matrix elements of the operators with dimension greater than five are not 
known. In order to justify retaining only the lowest dimension operators, the photon energy 
must be averaged over a region AE^ ^> Aqcd- At the present time these higher dimension 



operators introduce a significant uncertainty, since the photon spectrum is only measured 
experimentally over a region about 500 MeV from the end-point [|I|] . 

When operators in H e g other than O7 are included, the B — > X s 7 decay receives contri- 
butions from diagrams in which the photon couples to light quarks. It is well-known that for 
such processes, there are uncalculable contributions suppressed by a s , but not by powers of 
the scale associated with the process. Typically, the leading logarithms are calculable ||17|| , 



but terms suppressed by a logarithm (or equivalently by a s , but not by a power) can only 
be estimated using information on the fragmentation functions D q ^x and D g ^x deduced 
from other experiments or from models. While this may be worrisome, experience shows 
that usually the leading order perturbative QCD calculation provides an order of magnitude 
estimate of these effects ||18|| . Perturbative computations indicate that for weak radiative B 



decays into hard photons both the contribution of light quark loops ||, and the effects re- 
lated to decay functions of light partons into a photon || , are very smalLQ Therefore, these 
nonperturbative effects which are not power suppressed constitute less than five percent 
uncertainty in the theoretical prediction for the B — > X s 7 decay rate. 

There is no OPE that allows one to parametrize nonperturbative effects from the photon 
coupling to light quarks in terms of B meson matrix elements of local operators. Given 
this, it is perhaps not surprising that nonperturbative effects that come from the photon 
coupling to the charm quark contain B meson matrix elements of local operators that are 
suppressed by (Aqcd/^c) 2 rather than (AQ CD /mb) 2 . Recently, Voloshin identified such a 
nonperturbative correction to the B — > X s 7 decay rate ||15|| . This contribution arises from 



the interference of O2 with O-j corresponding to the diagram shown in Fig. 2, and can be 
studied using the operator product expansion. 



*For soft photons these effects are important. There are also interference effects where the photon 
couples to a light quark and to the charm quark, or to a light quark and through O7. These are 
also small for hard photons. 
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FIG. 2. Feynman diagram that gives rise to T 27 in Eq. (|j). Interchange of the photon and 
gluon couplings to the charm loop is understood. 

For a sufficiently heavy charm quark, nonperturbative corrections to the contribution of 
the interference of 2 and 7 to the decay rate can be computed from the discontinuity 
of the diagram in Fig. 2. Analogous diagrams with more gluons give effects suppressed by 
additional powers of Aqcd/^c- Denoting the gluon momentum by k, we work to all orders 
in k ■ q/m 2 c since the photon momentum q is of order m b , but neglect terms of order k ■ q/m b , 
k 2 /m 2 cb1 and m s /m b . The result of the loop integration is 



T 27 = --^-(B(v)\bm b a^q ^ J >f(l - 75 ) I, v b \B(v)} . 

ImB [m b v — q) z + it 

I^u is a complicated operator involving all powers of (q ■ iD)/m 2 . It is given by 
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(5) 



(2n + 4)! 

Here G Xv is the gluon field strength tensor and D denotes the covariant derivative. The 
contribution of T 27 to the B — > X s 7 decay rate is given by Eq. (||) with C 2 T 77 replaced by 
2C2C7 T 27 . 

For the leading n = term in Eq. (|J), the matrix element in Eq. (^) can be computed 
using the identity (19 



1 (B(v)\bT g s G aP b \B(v)) = ^ Tr {r (1 + f) a afi (1 + f)} 



2rriB 

valid for any Dirac structure V. The ratio of the decay rate from the n = term in Eq. 
to that from 7 is 

5T(B^X sl ) C 2 X 2 



(6) 



T(B^X sl ) 
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(7) 



The measured B* — B mass splitting gives A2 = 0.12 GeV 2 . Using this value for A2 and the 
values of C2 = 1.11 and C-j = —0.32 in Ref. 0, Eq. (|7|) implies that this O2 — O7 interference 
is about a three percent effect. This is an order of magnitude larger than the perturbative 
estimate of the contribution from the interference of O2 and O7 to the B — > X s 7 decay rate 
(which contains a gluon in the final state). 

The contribution of all terms in 1^ to Im T27 is 

ImT 27 = ^-t a " ( - 1 £< +3 r • • • r" (8) 

2m B ^ mj n+2 

x(B(v)\bT aP (q,v)(iD {lll ...iD lin} g a G a0 ) b \B(v)) S(q-v - 1/2) , 

where q = q/m b . T a/3 and a n are dimensionless and can be deduced from Eqs. (f|) and @. 
The indices fi\ . . . fi n are symmetrized, since they are dotted into g Ml . . . q^ n . Note that the 
derivatives iD^ . . . iD^ act on the gluon field G a p, and are determined by the spacetime 
dependence of the chromomagnetic field in the B meson. The n = term in Eq. (|8]) is a 
special case in that the (B(v) \ b T al3 g s G a p b \B(v)) matrix element is known from the B* — B 
mass splitting. The n = 1 matrix element vanishes by the equations of motion [F21]]. The 
n > 1 terms in Eq. (^) depend on an infinite series of unknown matrix elements. Estimating 

(B{v)\br#{q,v){iD fll ...iD ltn g a G a p)b\B{v))/{2m B ) ~ (Aqcd)" +2 , we see that the n > 1 
terms are "suppressed" compared to the n = term considered by Voloshin only by powers 
of m b A QCD /m 2 . 

In the limit where m c is fixed and rrif, — > 00, the higher order terms in Eq. (§) become 
successively more important and the expansion we have made is clearly inappropriate. (The 
whole sum in Eq. fl8|) is, up to logarithms, of order AqcD/^b-) in the limit where mf,/m c is 
held fixed and both masses become very large, the n > 1 terms in Eq. (|8|) are suppressed 
by powers of AQ CD /m c . Then the n = result, which is of order Aq CD /m 2 , dominates the 
sum. In the physical world, mbAq CD /m 2 ~ 0.6 (an equally reasonable estimate would be 
E™ x A/m 2 , which is also about 0.6). From Eq. (|5|) we see that ai/a = 4/15, g^/^o = 3/35, 
a 3/ a o = 16/525, etc., and asymptotically a n /a — > 3A/7r/(2™ +1 n 3//2 ) as n — > 00. The values of 
a n /a , n = 1,2,..., are small. This together with the asymptotic formula for large n suggests 



that the n > 1 terms in Eq. (|j) do not introduce a nonperturbative uncertainty greater than 
the value of the leading n = term. Nonperturbative effects from the interference of 0\ 
with Oj are expected to be smaller. 

Near the photon end-point region, another set of corrections become large. Expanding 
factors of iD that occur in the denominator of the strange quark propagator (these were 
neglected in Eq. (£D) yields corrections suppressed by powers of m*,. However, these correc- 
tions are proportional to derivatives of the delta function 8{E 1 — rrib/2), and they become 
as important as those in Eq. (||) in the end-point region. 

Consider next the contribution to the B — > X s 7 decay rate coming from the square of 
CiOi+C^C^- Diagrams like that in Fig. 2 should give a smaller nonperturbative contribution 
to the decay rate than the interference of O2 with O7 (i.e., these are order AQ CD /m 4 instead of 
order Aq CT) /m 2 ). But we know in this case that there is a contribution to the B — > X s 7 decay 
rate from B — > X s J/ip followed by J/ip — > 7 X, which is much larger than the perturbative 
calculation of the effect of (C\0\ + C2O2) 2 '. The combined branching ratio for this process 
is about 10~ 4 , while the perturbative estimate of the contribution of {C\Oi + C 2 2 ) 2 is less 
than 10 -5 . This might not present a serious difficulty for the comparison of experiment with 
theory, since the process B — > X s J/ip followed by J/ip -47I does not favor hard photons, 
and in any case it can be treated as a background and subtracted away. Further work on 
this issue is warranted. 

In this letter we examined uncertainties in the theoretical prediction for the weak ra- 
diative decay rate of B mesons into hard photons that come from nonperturbative strong 
interaction physics. We focused on effects that arise from photon couplings to light quarks 
and to charm quarks. For hard photons the first of these sources of theoretical uncertainty is 
less than five percent. This is smaller than the uncertainty in the Wilson coefficient C-?(mb) 
from uncalculated order a 2 s terms in its perturbative expansion. For the photon coupling 
to the charm quark, more work is needed to decide the size of the theoretical uncertainty 
associated with nonperturbative effects. 

The present experimental data on B — > X s 7 focuses on photon energies in the region 



Ej > 2.2 GeV 0. For comparison with this data, the largest theoretical uncertainty is 
from the contribution of higher dimension operators to the time ordered product T 77 which 
become more important in the end-point region. This uncertainty would be substantially 
smaller if the photon energy cut were reduced. 
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